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A recent paper by K. Xu et al. discussed a quan- 
tum depletion in gaseous Bose-Einstein condensation. 
In the ground state of an interacting Bose gas, not all 
the particles are in the lowest energy state, because the 
two-body interaction mixes the ground state components 
with atoms in the other states. For a homogeneous sys- 
tem the particle depletion is known well from Bogoliubov 
theory as 2] 770 = g^i/po^, where p is the atomic den- 
sity and a is the s-wave scattering length. The authors 
put the data of ^"^Na into above relation and say that "For 
^^Na at a typical density of lO^'^cm"^, the quantum de- 
pletion is 0.2%". Apparently, this is not a reliable value 
for the inhomogeneous system which a magneto-optical 
trap is applied. For an interacting and inhomogeneous 
system the difficulty comes from the fact that we do not 
know the every energy eigenvalue of the system except eo 
which is obtained by a macroscopic mean field wave func- 
tion so called Gross-Piataevskii equation(GPE). How- 
ever, one of the authors, W. Ketterle, already discussed 
the problem in an inhomogeneous system seriously and 
obtained an effective method for that We'll apply 
this method to obatin a better prediction. 

We first consider a 3D system of N non-interacting 
Bosons confined by a harmonic potential of angular fre- 
quency oj. When the Bose distribution function in the 
grand canonical ensemble at temperature T is /^(e) and 
the density of states of the isotropic 3D harmonic oscil- 
lator is p{e) , the number of particles at temperature T is 
obtained by / fB{e)p{e)de 0,11,1111. It is 

N^N,+ ('^)\,i^) + l('^)\.i.) + ..., (1) 
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where gn{z) = Y^^i -^V" is the Bose functions, z = 
g/3(/i-£o) jg tj^g fugacity between and 1, Eq = (3/2)?iu; 
is the zero-point energy of the harmonic oscillator, and 
/i(T) is the chemical potential. Aq is the ground state 
occupation number given by Aq — l/{z~^ — 1). The 
others are common notations. 

For interacting system there is a shift of the zero point 
energy by the interaction. Then, the new ground state 
occupation number is written as 

(2) 
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where U — U{a,No) is the interaction energy between 
atoms and may obtained from the solution of the GPE. 
We consider the positive scattering length only, so does 


U > 0. Then, the condensate fraction of the interacting 
system is definitely lower than that of the ideal system. 

The chemical potential p{T) is obtained numerically 
by the following way. The transition temperature Tc is 
determined that at the onset of condensation A'o 
and z 1 from Eq. (Q. It is 
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where Tq = {Tiuj/kB) {A^/C(3)} and C, is the Riemann- 
zeta function. = 328 - 713nA: for iV = 10^ ~ 10^ 
and uu = lO'^sec^^. The condensate fraction is obtained 
from Eq. as a function of temperature 


Ao^^ /T 
N \To 
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The author's prediction of the quantum depletion 0.2% 
provides that the transition temperature should be below 
40 - 90nK for A^ = 10^ - 10^. Actually, the author's 
estimated value 770 = 0.2% is too small for the inhomo- 
geneous system to be accepted. We obtain the fugacity 
z{T) first by solving Eq. numerically, and then ob- 
tain the chemical potential p{T). We use this p{T) for 
the estimation of A"™*. 

Comparing A^o and N™*, the over all behavior of the 
condensate fraction in Eq. Q as a function of scattering 
length and temperature can be estimated by the approx- 
imation 


Af(5"*(a,T) - A^o(T)e-^'^(°'^°l 


(5) 


The ground state occupation number falls exponentially 
as the interaction increases. 
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